Spontaneously generated X-shaped light bullets 
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We observe the formation of an intense optical wavepacket fully localized in all dimensions, i.e. 
both longitudinally (in time) and in the transverse plane, with an extension of a few tens of fsec and 
microns, respectively. Our measurements show that the self-trapped wave is a X-shaped light bullet 
spontaneously generated from a standard laser wavepacket via the nonlinear material response (i.e., 
second- harmonic generation), which extend the soliton concept to a new realm, where the main 
hump coexists with conical tails which reflect the symmetry of linear dispersion relationship. 
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Defeating the natural spreading of a wavepacket (WP) 
is a universal and challenging task in any physical context 
involving wave propagation. Ideal particle-like behavior 
of WPs is demanded in applications, such as microscopy, 
tomography, laser-induced particle acceleration, ultra- 
sound medical diagnostics, Bose-Einstein condensation, 
volume optical-data storage, optical interconnects, and 
those encompassing long-distance or high-resolution sig- 
nal transmission. The quest for light WPs that are both 
invariant (upon propagation) and sufficently localized in 
all dimensions (3D, i.e., both transversally and longitu- 
dinally or in time) against spreading " forces" exerted by 
diffraction and material group- velocity dispersion (GVD, 
k" = d 2 k/duj 2 \ UJo ) has motivated long-standing stud- 
ies, which have followed different strategies in the linear 
[J E>] and nonlinear 0, Q regime, respectively. 

In the linear case, to counteract material (intrinsic) 
GVD, one can exploit the angular dispersion (i.e., de- 
pendence of propagation angle on frequency) that stems 
from a proper WP shape. The prototype of such WPs is 
the X-wave Q, a non- monochromatic, yet non-dispersive, 
superposition of non-diffracting cylindrically symmetric 
Bessel Jo (so-called conical or Durnin [l|) beams, ex- 
perimentally tested in acoustics Q, optics Q and mi- 
crowave antennae [5| . Importantly, in the relevant case of 
WPs with relatively narrow spectral content both tempo- 
rally (around carrier frequency ujo) and spatially (around 
propagation direction z, i.e. paraxial WPs), X- waves 
require normally dispersive media (k" > 0). In this 
case, a WP with disturbance E(r, t, z) exp(ikoz — iuo^T) 
(ko = k(uo), r 2 = x 2 -f ?/ 2 ), has a slowly- varying envelope 
E = E(r, t, z) obeying the standard wave equation 
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Laplacian, where = d 2 r + r~ x d r is the transverse 
Laplacian, and we limit our attention to luminal WPs 
traveling at light group- velocity 1/k' = dk/duo\~^ by 
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FIG. 1: Features of normally dispersive media (k" > 0): (a) 
spatio-temporal dispersion relation associated with paraxial 
wave equation JJJ; Different representations [(b) sections of 
color isointensity surfaces \E(x, y, t)\ 2 =const.; (c) intensity 
vs. r, t] of a non-diffractive, non- dispersive linear X-wave 
(A = 10 fs, k" = 0.02 ps 2 /m, k = 10 7 m" 1 ). 



introducing the retarded time t = T — k'z in the 
WP bary centre frame. Propagation-invariant waves 
E(r, t, z) = E(r, t, z = 0) exp(i/3z) can be achieved when- 
ever their input spatio-temporal spectra E(K, Vt,z = 0) 
lie along the characteristics of the dispersion relationship 
k"n 2 /2 - K 2 /(2k ) = (3, which follows from Eq. Q in 
Fourier space {K, Q) (K is the transverse wavevector re- 
lated to cone angle with z-axis c± sinf9 = K/ko, and 
Q = co — (Jo)- In the normal GVD regime (k" > 0) 
these curves, displayed in Fig. ^a), reflect the hyper- 
bolic nature of the wave equation JQ and show the com- 
mon asymptotic spectral X-shape associated with the 
lines K = ±2kok"fl ((3 = 0). When the spectral com- 
ponents lying along such X are superimposed coherently 
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(in phase), the field in the physical space (r,i) also re- 
tains a propagation invariant X-shape, as entailed e.g. 
by the exact solution E = Re{[(A - it) 2 + /coZc'V 2 ]" 1 / 2 } 
of Eq. JU shown in Fig.^b-c). Here A represents dura- 
tion of the X-wave central hump. Main features of these 
waves are the conical (clepsydra) 3D structure and the 
slow spatial decay (1/r characteristic of Jo components 
1]), displayed by Fig.^b) and (c), respectively. 

Conversely, in the nonlinear (high intensity) regime, 
non-spreading WPs exploit the idea that, in self- focusing 
media, the nonlinear wavefront curvature can simulta- 
neously balance the curvature due to diffraction and 
to GVD, combining features of spatial |£j and tempo- 
ral [HI solitons to form a bell-shaped 3D-localized WP 
E(r,i), so-called light-bullet In sharp contrast with 
X-waves, such compensation strictly requires anomalous 
GVD (k" < 0) 6], thus implying that along the WP tails, 
where Eq. (1) still holds true, space r and time t play the 
same role giving rise to strong WP localization pj| . Sta- 
ble trapping, however, has been observed only in setting 
of reduced dimensionality (2D) including also other 

contexts, e.g. spin [ll[ or atomic waves (3D kinetic 
energy and atom-atom attractive interactions act exactly 
as diffraction- GVD and self-focusing, respectively) where 
similar trapping mechanisms hold true. 

In this work, we outclass the two approaches by demon- 
strating that space-time localization in the normal GVD 
regime becomes accessible in the nonlinear regime. Trap- 
ping is accomplished by mutual balance of intrinsic, 
shape-induced, and nonlinear contributions in a new type 
of WP, namely a nonlinear X-wave, which permits to get 
over two limitations at once. First, in contrast with lin- 
ear X-waves, whose observation requires non-trivial input 
beam shaping 0, 0, our experiment reveals a remark- 
able " mode- locking" process that, starting from a con- 
ventional (gaussian) laser WP, spontaneously performs 
the reshaping into a localized X-shaped WP. Second, we 
believe this to be the first genuine nonlinear trapping in 
full-dimensional 3D physical space, since to date mate- 
rial and/or instability limitations |y,L| have rendered the 
observation of light 3D bullets elusive. 

Figure El describes the strong localization features ob- 
served after propagation in a 22 mm long sample of 
lithium triborate (LBO) crystal. At the input 

we launch a laser WP at fundamental frequency (FF) 
ujq = 27rc/Ao, Ao = 1060 nm, with gaussian profile in 
both t (with FWHM duration in the 100 - 200 fs range ) 
and r (45 /im FWHM at waist, located few mm before the 
crystal so that the input beam is slightly diverging) . The 
LBO crystal is tuned for generation of optical second- 
harmonic (SH) in the regime of relatively large positive 
phase mismatch Afc = 2k(wo) — k(2coo) = 30 cm -1 or 
effective self- focusing for the FF beam [IBj. When the 
input energy exceeds about 0.25 /iJ, mutually trapped 
localized WPs at FF and SH are observed (see Fig. EJ). 
Time-integrated measurements of spatial profiles (Fig. [21 
top) indicate that diffraction is fully defeated to yield 
a spatial soliton-like beam 0, while temporal autocor- 



FIG. 2: Top: output transverse beam profiles in the linear 
(white) and nonlinear (green: SH; red: FF) regime. The 
diffracted (left) and localized (right) spots recorded at the 
output of the crystal are also shown. Bottom: temporal auto- 
correlation traces (collinear technique) showing the transition 
from the linear regime (red) to the compressed WP (green). 

relations reveal a single pulse which is strongly com- 
pressed (down to ~ 20 fs, see Fig. bottom). Re- 
markably nonlinear mixing balances the two highest- 
order dispersive effects, namely the tendency of FF and 
SH to walk-off due to group-velocity mismatch (GVM) 
SV = k' '(2coo) — k'(uJo) an d GVD, which act in the linear 
regime over characteristic length scales 50 and 2.5 times 
shorter than the crystal length, respectively. Relying on 
nonlinear scenarios known to date, the present result is 
unexpected. In fact, the normal GVD and the strong 
GVM of LBO do not allow an explanation in terms of 
light bullets, whereas the dynamics of self- focusing with 
normal GVD is dominated by pulse splitting without en- 
visaging localization whatsoever |l4j . 

In order to get a deeper understanding, we performed 
a set of numerical experiments by integrating the well- 
known model 15], which generalizes Eq. to nonlinear 
coupling of envelopes E m (r, t, z) (m = 1,2, carrier ttiuoq) 

L(2^ ) E 2 + SVd t E 2 + x&i exp(zAZcz) = 0, [ } 

In Fig. (top) we show snapshots of the space-time evo- 
lution. It is clear that, in the first stage (z = 10 — 15 mm), 
the generation of light at SH is accompanied by strong 
self-focusing, in turn inducing pulse compression (in spite 
of the normal GVD that would cause temporal broaden- 
ing of plane- waves). The WP becomes asymmetric since 
self-focusing is stronger where the SH, which lags behind 
because of GVM, tends to be. More importantly, during 
this process the WP undergoes a strong reshaping, where 
a large fraction of energy is radiated off-axis to form the 
tails of a conical wave (see snapshot at z = 20 mm). Af- 
ter this stage, however, the collapse stops and the WP 
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propagates (locked with the SH) with immutable shape 
and nearly constant energy, duration, and size. 

By changing the parameters in Eqs. we can con- 
clude that: (i) nonlinearity is the key element that drives 
the reshaping and holds the WP together. In fact, by 
switching it off after the transient (i.e., x = f° r z > 20 
mm), the WP exhibits strong diffraction and extremely 
fast FF-SH walk-off. The reshaped WP can by no means 
be considered a linear X-wave; (ii) in our crystal, GVM 
is the dominant dispersion term affecting the asymptotic 
duration and width of the localized WP (both decrease 
for smaller GVM). However, (symmetric) X-shaped WPs 
are formed also in the ideal case of vanishing GVM, 
provided that GVD is normal; (iii) the WP reshaping 
strongly affects the phase modulation process, which re- 
veals the contribution of an effective anomalous GVD; 
(iv) by including additional cubic nonlinearities 0, the 
phenomenon remains qualitatively unchanged. 

A more rigorous ground for explaining the dramatic 
reshaping shown in Fig. |3| is offered by the investiga- 
tion of two different aspects, both implicitly accounted 
for by Eqs. J2J): (i) the stability of continuous plane- 
waves; (ii) the existence of nonlinear X-shaped eigen- 
solutions. Though somewhat idealized, the results of 
this analysis allow gathering the two basic features of 
the spatio-temporal dynamics, namely the transient re- 
shaping of the input gaussian-like WP followed by the 
quasi- stationary regime. 

The linear stability analysis of z-invariant solutions of 
Eqs. J2J) with ideally vanishing spatio-temporal spectral 
width (i.e., continuous plane waves) always reveal the 
presence of exponentially growing weak (up to noise level) 
perturbations with definite frequency K, Q. However, 
the instability features reflect the symmetry of the linear 
wave equation, thus being qualitatively different in the 
normal and anomalous GVD regime, respectively. While 
anomalous GVD leads to narrow bandwidth features as 
in conventional spatial or temporal modulational insta- 
bility 0, the hyperbolic structure of the diffraction- 
dispersion operator in the normal GVD regime of our 
experiment leads to exponential amplification of conical 
wave perturbations (Bessel Jo beams) with frequencies 
basically approaching the asymptotes of Fig. [IJa) [Ti| . 
When growing spontaneously from noise, the amplified 
components in the virtual infinite bandwidth (in real- 
ity limited by non-paraxiality and higher-order dispersion 
not accounted for in Eqs. (1-2)) leads to colored conical 
emission. Our calculations show that the phenomenon 
persists under dynamical conditions (unseeded SH gen- 
eration) and when pumped by short-pulse narrow-beam 
inputs, as in our experiment. In the latter case, the am- 
plification of proper frequency components of the WP 
preserves the mutual phase coherence, thus acting as a 
trigger which drives the transformation of the WP into 
the X-wave shown in Fig. [31 In other words, it is the 
conical instability that probes the symmetry of the un- 
derlying linear system in amplifying those components 
which allow both diffraction and dispersion of the whole 
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FIG. 3: Numerical simulation of the localization process ob- 
served in our LBO sample [from Eqs. with k" '(ooo) = 
0.016 ps 2 /m, k"(2u ) = 0.089 ps 2 /m, SV = 45 ps/m, 
X = 7x 10" 5 W" 1/2 (deff = 0.85 pm/V)] from an input 
gaussian WP with 45 /xm waist, 170 fs duration, and 0.26 /iJ 
energy. Top: snapshots showing the FF beam in the t — x 
plane. Bottom: evolution of the beam diameter, duration, 
and energy content of main localized hump. The dashed line 
refers to more energetic (e.g., slightly wider and longer) in- 
put WPs for which even the small residual diffraction (solid 
line) is removed. In this case, however, the reshaping leads to 
higher peak intensity close to damage threshold of LBO. 



WP to be removed. 

A second, strong argument in favour of the nearly 
asymptotic character of the evolution shown in Figs. 03 
is the existence of stationary localized solutions of 
Eqs. (0). We have recently shown, indeed, that 
the natural propagation-invariant [i.e., with dependence 
E m (r, t) exp(im(3z), m = 1,2] localized eigensolutions of 
SH generation process in the normal GVD regime (either 
with or without GVM) are indeed nonlinear X waves ^3 • 
These waves are similar to those shown in Fig. E^b-c) 
except for the fact that their peak intensity is related 
to their duration and width through the nonlinearity. 

Since the measurements in Fig.[21give information only 
about the WP central hump, in order to confirm experi- 
mentally that the observed strong focusing and compres- 
sion dynamics is indeed driven by the formation of an 
X-wave, we have performed different additional measure- 
ments. First, we have characterised the propagation of 
the WP in air after the LBO nonlinear crystal. We ob- 
serve that the beam diffracts less than a gaussian beam 
of the same width (the divergence angle is 2.5 times less), 
and this net sub-gaussian diffraction witnesses a Bessel- 
like feature of the WP. Moreover, temporal broadening 
to 150 fs after only 10 cm of propagation in air is an indi- 
cation that the WP develops strong angular dispersion. 
Both observations are in good quantitative agreement 
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FIG. 4: (a) Output spatio-temporal intensity profile, as mea- 
sured in air, 5 mm from the crystal output face. Insets: trans- 
verse intensity pattern in (x, y) plane measured at peak (t = 
fs) and with t = 70 fs delay, (b) As in (a), numerical result 
from Eqs. (2). (c) As in (b) calculated right on crystal output. 

with data obtained numerically from Eqs. (0). However, 
it is only the tomography of the output WP, i.e. map- 
ping the WP intensity in space and time, that can give 
direct unequivocal evidence for the formation of an X- 
wave. To this end we have developed a new technique 
based on an ultrafast nonlinear gating, or a scanning 
cross-correlation technique, realized by frequency mixing 
the WP under investigation with a 20 fs, high contrast, 
steep front, probe pulse in the visible, which is uniform 
over a large (few mm 2 ) area (the details of the set up 
will be presented elsewhere). Thanks to the high quality 
of the probe and the use of a very thin (20 fim) BBO 
mixing crystal, the apparatus has high temporal resolu- 
tion. The intensity map of the output WP is reported 
in Fig. Ufa). The insets also show the measured beam 
profile in the transverse plane at time t = (peak) and 
t = 70 fs (far from peak), respectively. The measured 
profile clearly shows the features of an X-wave with a 



conical structure, which emanates from a strongly local- 
ized central spatio-temporal hump. Unlike conventional 
pulse splitting ^4)> nere splitting occurs only sufficiently 
off-axis (x ~ 100/im) in the WP low-intensity portion. 
For comparison we also report in Fig. Efb) the profile 
calculated from Eqs. (2) under the same conditions, the 
agreement being excellent. The fringe-like structure that 
appears for large delays in both Fig. 0|a) and (b) is due 
to 5 mm of free-space propagation in air outside the crys- 
tal. Although the calculated profile on the output face 
of the crystal indeed shows that such fringes disappear 
[see Fig.0Jc)], measurement with perfect imaging on the 
output LBO face reveals saturation of the mixing process 
due to a too intense peak. 

In summary, we have reported the first evidence that 
the natural 3D (temporal and spatial) spreading of a fo- 
cused ultrashort wavepacket can be balanced in trans- 
parent materials at high intensity. The underlying mech- 
anism is the spontaneous formation of a X-wave charac- 
terized by an intense (i.e., nonlinear) central hump self- 
trapped through mutual balance with (essentially linear) 
dispersive contributions associated with coexisting slowly 
decaying conical tails. While our experiment is carried 
out by exploiting self-focusing nonlinearities arising from 
quadratic nonlinearity, we envisage the general role that 
self-trapping mediated by nonlinear X- waves can have for 
a wide class of materials and applications encompassing 
centrosymmetric optical (Kerr) media Q, Bose-Einstein 
condensation and acoustics 0. 
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